Evidence for site-centered stripes from magnetic excitations in CuO superconductors 
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The success of models of coupled two-leg spin ladders in describing the magnetic excitation spectrum of 
La2_vBaj-Cu04 has been widely interpreted as evidence for bond-centered stripes. Here, we determine the 
magnetic coupling induced by the charge stripes between bond- or site-centered spin stripes modeled by two- 
or three-leg ladders, respectively. We find that only the site-centered models order. We further report excellent 
agreement of a fully consistent analysis of coupled three-leg ladders using a spin wave theory of bond operators 
with the experiment. 

PACS numbers: 74.72.-h, 74.20.Mn, 75.10.-b, 75.25.+Z 



The mechanism of superconductivity in the presence of lo- 
cal antiferromagnetism in copper oxides is considered an out- 
standing problem in contemporary physics''^. The materials 
are described by mobile charge carriers (holes) doped into 
a quasi-twodimensional spin 1/2 antiferromagnet^. Inelas- 
tic neutron scattering experiments have revealed a magnetic 
resonance peak'*'^ and, in some compounds, periodic mod- 
ulations in the spin and charge density (stripes)'^'^. Tran- 
quada et al}^ found that the magnetic excitation spectrum 
of stripe ordered Lai 875Bao.i25Cu04 looks similar to disor- 
dered YBaiCuaOfi+.i^ or Bi2Sr2CaCu208+^'"*, and observed 
that the data are consistent with bond-centered stripes mod- 
eled by two-leg ladders. This experiment, and its interpre- 
tation, is considered of key importance for the field of high- 
temperature superconductivity. The reason these findings are 
so vigorously discussed is that they may provide the decisive 
hint as to within which framework these conceptionally sim- 
ple yet inscrutable systems are to be understood, and hence 
ultimately lead to a complete theory. 

At present, there is no consensus with regard to such a 
framework, but instead a fierce competition among different 
schools of thought. One of these schools'^ '^'" '^ attributes the 
unusual properties of the doped, two-dimensional antiferro- 
magnets to their propensity to form stripes, or their proximity 
to a quantum critical point (QCP) at which stripe order sets 
in. The resulting picture is highly appealing. Static stripes 
have been observed^- only in certain compounds, most no- 
tably La2-ASrvCu04 at a hole doping concentration x - ^, and 
are known to suppress superconductivity. On the other hand, 
the mere existence of stripes would impose an effective one- 
dimensionality, and hence provide a framework to formulate 
fractionally quantized excitations. This one-dimensionality 
would be roughly consistent with an enormous body of exper- 
imental data on the cuprates, including the electron spectral 
functions seen in angle-resolved photo emission spectroscopy 
(ARPES). The charge carriers, the holons, would predomi- 
nantly reside in the charge stripes, as they could maximize 
their kinetic energy in these antiferromagnetically disordered 
regions. In the spin stripes, by contrast, the antiferromagnetic 
exchange energy between the spins would be maximized, at 
the price of infringing on the mobility of the charge carri- 
ers. Most importantly, the spin stripes would impose a cou- 
pling between the charge stripes, which would yield an effec- 
tive, pairwise confinement between the low-energy spinon and 



holon excitations residing predominantly in the charge stripes. 
The mechanism of confinement would be similar to that of 
coupled spin chains or spin ladders'^ The holes would 
be described by spinon-holon bound states, and the domi- 
nant contribution to the magnetic response measured in Tran- 
quada's as well as all other neutron scattering experiments 
would come from spinon-spinon bound states. 

The similarity of the "hour-glass" spectrum shown in 
Fig. 4b of Tranquada et al}^ (which is reproduced for com- 
parison in Fig.[T}3) and the "elephants trousers" observed by 
Bourges et al.^ (see Fig. 3 of their manuscript) provides the 
most striking evidence in favor of the picture advocated by this 
school, which attributes the anomalous properties of generic, 
disordered CuO superconductors to the formation of dynamic 
(rather than static) stripes, which fluctuate on time scales slow 
compared to the energy scales of most experimental probes. 
This picture received additional support by Xu et al}^, who 
observed that the magnetic response of Lai g75Ba() 125CUO4 
at higher energies is independent of temperature, while the 
stripe order melts at about T^i ~ 54K. Measurements on 
'untwinned' samples of YBa2Cu306.6, where one would ex- 
pect the dynamical stripes to orient themselves along one of 
the axis, however, exhibit a strong anisotropy in the response 
only at energies below the resonance'^, while the response 
is fourfold rotationally symmetric at higher energies^*^. This 
may indicate that the formation of stripe correlations, be it 
static or dynamic, is a low energy phenomenon, while the high 
energy response probes itinerant antiferromagnets at length 
scales on which the stripes are essentially invisible. Another 
extremely appealing feature of the experiment by Tranquada 
et al}^ is that it immediately suggests a model of ferromag- 
netically coupled two-leg ladders, as the upper part of the 
measured spectrum agrees strikingly well with the triplon (or 
spinon-spinon bound state) mode of isolated two-leg ladders 
(see Fig. [TJj). The experiment hence appears to point bond- 
centered rather than site-centered stripes {i.e., stripes as de- 
picted in Fig.|2K rather than Fig.|2jJ), and thereby to resolve 
a long outstanding issue within the field of static stripes. 

This interpretation received support by theoretical stud- 
ies^'"^''. Vojta et al.^^ used a bond operator formahsm^i to 
study a spin-only model of coupled two-leg ladders (as de- 
picted in Fig. |2l\) with J\\ - J, took into account a bond- 
boson renormalization of J, and assumed a value J' for the 
ferromagnetic coupling between the ladders which is large 




FIG. 1: (Color online) Comparison of dispersions resulting from dif- 
ferent theories with experimental data, a, Superpositions of cuts 
along (fcj, Ti) and {n, k,j) for the lowest mode of the SWT with J = 
140 meV and J' = 0.07 J for site-centered stripes described in the 
text (red) superimposed with the experimental data obtained by in- 
elastic neutron scattering'^ (black), b, The neutron data, with a 
triplon dispersion of a two-leg ladder superimposed (red line), re- 
produced from Tranquada et alM (Reprinted by permission from 
Macmillian Publishers Ltd: Nature 429: 534-538, ©2005). 



enough to close the spin gap of the ladders, i.e., to induce long 
range magnetic order. Within their approximations, a value of 
J' - -0.067 is sufficient. This value is not consistent with 
previous studies^^"^^, but as no method to calculate or even 
estimate the true J' induced by charge stripes was available, 
it did not seem a problem at the time. The spectrum they ob- 
tained agrees well with experimental data measured by Tran- 
quada et alM, and hence appeared to justify their assump- 
tions a posteriori. They concluded in favor of bond-centered 
stripes. This conclusion was independently strengthened by 
Uhrig et al.^-, who used the method of continuous unitary 
transformations to study a model of ferromagnetically cou- 
pled two-leg ladders, and observed that the critical value of 
J'^ can be significantly reduced if a cyclic exchange term 7cyc 
on the ladders is included^^. They likewise fine-tuned /' to 
the QCP where the gap closes and long-range magnetic order 
ensues, and reported good agreement with experiment. Sei- 
bold et fl/.-^'^" calculated the magnetic response for a range 
of dopings within the time-dependent Gutzwiller approxima- 
tion, and found good agreement with the measured data for 
both bond- and site-centered stripe models. 

The first question we wish to address here is whether the 
key assumption that J' is large enough to induce order is valid. 
There exist several estimates for the critical value J[. required 
if the coupling between isotropic ladders is antiferromagnetic. 
Gopalan et ali^ find J[. x 0.257 in a simple mean-field treat- 
ment of bond-bosons. Quantum Monte Carlo (QMC) calcula- 
tions by Tworzydlo et al.^^ yield J'^ - 0.30(2)7, a value sub- 
sequently confirmed by Dalosto et al.^''. We have redone the 
mean-field calculation of Gopalan et al.^^ for ferromagnetic 
(FM) couplings J'^ < 0, and find that within this approxima- 
tion, the absolute value of J'^ is independent of the sign of the 
coupling. QMC calculations by Dalosto et al.^^, however, in- 
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FIG. 2: Models of spin stripes: a. Bond-centered stripes modeled 
by two-leg ladders and b. Site-centered stripes modeled by three- 
leg ladders, c. Definition of bosonic creation operators in the eight- 
dimensional Hilbert space spanned by a rung on sublattice of the 
ladders shown in b. 



dicate that the true value is at least = -0.47 (see Fig. 6b 
of their article). The physical reason why a significant cou- 
pling between the ladders is required to induce magnetic order 
is that the individual two-leg ladders possess a gap of order 
A Si 7/2. As a cyclic exchange term 7cyc ~ 0.257 reduces 
this gap by a factor of two2i, we would expect that 7^ would 
likewise be reduced by a factor of two. We hence conclude 
that a FM coupling of at least somewhere between 7^ = -0.27 
and -0.47 is required, depending on the strength of a possible 
cyclic exchange term. 

The explicit calculation of the ferromagnetic coupling in- 
duced by the charge stripes between the spin stripes we de- 
scribe below, however, yields 7' - -0.057. The coupling is 
hence insufficient to induce order in a model of bond centered 
stripes. 

For a model of site-centered stripes described by antiferro- 
magnetically coupled three-leg ladders, as shown in Fig. |2jJ, 
the critical coupling required for long range ordering to set 
in is by contrast 7^ = 0. The reason is simply that there 
is no need to close a gap, as the three-leg ladders are in- 
dividually gaplessi^. A conventional spin wave analysis for 
such a spin-only model of three- and four-leg ladders was per- 
formed by Yao et al.^\ who found that their approximation 
agrees reasonable well with the experimental data if they take 
7' - 0.057 and 7' = -0.097 for coupled three- and four- 
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FIG. 3: Finite size geometries with unfrustrated periodic boundaries 
for a, bond- and b, site-centered stripe models. The spin stripes are 
localized by a staggered magnetic field B as indicated by the signs in 
the grey shaded areas. 



leg ladders, respectively. The calculation we present in the 
following, however, singles out J' - 0.07 J for the antiferro- 
magnetic coupling between spin stripes modeled by three leg 
ladders. 

To determine the effective coupling J' between the two- and 
three-leg ladders representing bond- or site-centered stripes, 
respectively, we have exactly diagonalized 16 site clusters of 
itinerant spin 1/2 antiferromagnets described by the (nearest- 
neighbor) t-J modeP with J - 0.4?, two holes, and peri- 
odic boundary conditions (PBCs), in which the stripes are 
locaUzed through a staggered magnetic field B in the grey 
shaded areas shown in Figs.|3}\ and|3l}. We then compare the 
ground state energies we obtain for clusters with the unfrus- 
trated PBCs shown in Fig.[3]with the ground state energies we 
obtain for clusters with frustrated PBCs, in which the 16-site 
unit cells shown on the right are shifted by one lattice spacing 
to the top, such that sites 15 and 1, 16 and 2, etc. are nearest 
neighbors. We then consider spin-only Heisenberg models of 
two- and three-leg ladders (consisting of only the sites in the 
shaded areas in Figs.[3]\ and[3j3) subject to the same staggered 
field B and couple them ferromagnetically or antifeiTomagnet- 
ically by J', respectively, as indicated. We again compare the 
ground state energies for unfrustrated PBCs, where J' cou- 
ples sites 6 and 1, 7 and 2, etc. for the two leg ladders shown 
in Fig.[3K, with frustrated PBCs, where J' couples sites 7 and 
1, 8 and 2, etc. Finally, we determine J' such that the dif- 
ference in the ground state energies between frustrated and 
unfrustrated PBCs in the t-J clusters matches this diff'erence 
in the spin-only ladder models. 

With B = 0.2257 and B = 0.1707 for the bond- or site- 
centered stripe models we obtain J' = -0.0517 and 7' - 
0.01 IJ, respectively. The values for B are chosen such that the 
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FIG. 4: (Color online) Constant energy slices of the neutron scat- 
tering intensity for the lowest mode of the SWT with 
J = 140 meV and 7' = 0.07 7. We have replaced the (5-functions in 
frequency by Lorentzians with half-width A = 0.05 7 and averaged 
over both stripe orientations (i.e., horizontal and vertical). 



magnetic localization energy E^ag. - -B'Zji(~^yS~i is equal 
for both types of stripes, and such that the mean value of the 
staggered magnetization {-S\ -H 5^) in the t-J cluster for the 
site-centered stripe matches the value we obtain in the spin 
wave theory (SWT) of three-leg ladders described below. We 
estimate the eiTor resulting from finite size coiTections and the 
use of a nearest-neighbor t-J model to be of order ±0.017. 

In the remainder of this Rapid Communication, we will 
show, albeit not in full detail, that a fully consistent SWT 
of bond operators representing the eight-dimensional Hilbert 
spaces on each rung of the three-leg ladders agrees perfectly 
with the experimental data if and only if the coiTect, calcu- 
lated value 7' = 0.077 is used for the coupling. This discus- 
sion will complete our argument showing that the experiment 
by Tranquada et al. provides evidence for site-centered, and 
not bond-centered, stripes. 

To begin with, consider a single rung on sublattice J[ with 
three spins S\,S2, and S3, as indicated in Fig. |2j3. Diagonal- 
ization of the Hamiltonian - J{s\S2 + $2*3) yields a spin 
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doublet {p-1/2) , pi/2)} with energy E - -J and a quadruplet 

c-3/2) , |c-i/2) , |ci/2) , 1^3/2)} with E - J, both of which are 
symmetric under mirror reflection (interchange of sites 1 and 
3), and another doublet {|fl-i/2) , |fli/2)} with E — Q, which is 
antisymmetric under this mirror reflection. The subscripts in 
the kets label the eigenvalues of s[^^. We define a fiducial state 
1^-1/2) for sublatticeyi via |&_i/2) = |/7_i/2) cos 0+|c_i/2) sin 
which interpolates between the s\^^ - —j quantum ground 
state l/j-1/2) of the isolated rung for <p = and the classically 
Neel ordered state |iti> for (p - arctan(-^) = 0.6155. We 

then introduce bosonic operators = |a-i/2)(S-i/2l ^tc., as 
indicated in Fig.[2t. Written in terms of these operators, the 
Hamiltonian //'^ of a single rung will contain a linear term 
|(cj + CQ)sin2(^. Finally, we introduce a similar formalism 
with operators Aq = |fli/2){^i/2l etc. and a fiducial state IS1/2) 
with 5;^, = +5 for the rungs on sublattice S, and express all 
the spin operators on all the individual sites in terms of the 
bosonic operators. 

When we couple the rungs with intra-ladder couplings - 
J and inter-ladder couplings J' (see Fig. |2j3), we keep terms 
up to quadratic order in the bosonic operators, and adjust the 
angle (p introduced above such that the terms linear in the op- 
erators vanish. This yields (p = 0.341 for J' = 0.077. We in- 
troduce Fourier transforms of our bosonic operators into mo- 
mentum space using the unit cell indicated by the grey area in 
Fig. |2jJ, and rewrite the total Hamiltonian in terms of those. 
As the linear SWT preserves both the mirror reflection sym- 
metry and the total spin quantum number 5 -', the Hamiltonian 
separates into several terms, H - E^ + Hao + HfO + ^01 + 
Hc2 + Hbi,ci,c-i, where £0 is a contribution to the ground state 
energy, Hao contains only the operators a^, Aq, Aq, and A,j, 
and so on. The low-energy spectrum we are interested in is 
exclusively contained in 

i/M,.u-i=X(*I^A-5tr(i/4 



where = b\ ,, C_^ ,, c^^ ,, c_^ _f and is a 

6x6 matrix. We diagonalize Hhi ,-\,c-i via a six-dimensional 
Bogoliubov transformation at each point in k space. 

Cuts of the dispersion of the lowest mode we find are 
shown in Fig. [T]\. They agree strikingly well with the super- 
imposed experimental data by Tranquada et al.—. Constant 
energy slices of the neutron scattering intensity ;(f^"(^,w) of 
this mode, which may be compared to the data reported in 
Fig. 2 of Tranquada et ali^, are shown in Fig. |4] The other 
two modes of Hi,\^c\,c-\^ as well as all modes from the other 
terms in ft, are only weakly dispersing and start at energies of 
order 27, i.e., in a regime where we would no longer expect 
our spin wave theory to yield reliable results. 

Finally, let us remark that for 7' < 0.57, the saddle point 
energy in our model is given to a highly accurate approxima- 
tion by ojin, n) K \ .47 V7'7. Since we expect 7' to decrease 
with decreasing doping, we would expect a similar doping de- 
pendence for the saddle point energy.22^. 

In conclusion, we have shown that models of coupled two- 
leg ladders describing bond-centered stripes cannot explain 
the magnetic spectrum of La2-ABa.tCu04'^ as the coupling 
induced by the charge stripes between the ladders is insuffi- 
cient to induce long-range magnetic order We have further 
shown that a model of coupled three-leg ladders describing 
site-centered stripes accounts accurately for the experimen- 
tal data. The experiment hence provides evidence for site- 
centered, and not bond-centered, stripes. 

We wish to thank M. Vojta, T. Ulbricht and P. Wolfle for nu- 
merous discussions. This work was supported by the German 
Research Foundation under grant FOR 960. 
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